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Abstract. We report new results on the non-conserved dynamics of paraUel steps on vicinal surfaces in the 
case of sublimation with electromigration and step-step interactions. The derived equations are valid in the 
quasistatic approximation and in the limit /^^ ^ Id ^ l± ^ h, where / is the inverse electromigration 
length, Id the diffusion length, l± the kinetic lengths and h the terrace widths. The coupling between 
crystal sublimation and step-step interactions induces non-linear, non-conservative terms in the equations 
of motion. Depending on the initial conditions, this leads to interrupted coarsening, anticoarsening of step 
bunches or periodic switching between step trains of different numbers of bunches. 



PACS. XX.XX.XX No FAGS code given 
1 Introduction change of the crystal volume is a result of the movement 
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f^y^ of the steps. On this scale we can reduce a surface with 

For the theoretical study of homoepitaxial growth and 

. . . . . straight steps to a one-dimensional step train. Such a sur- 

I sublimation of a crystal in contact with the gas phase it 

>' ' ' . . . . . . face may undergo step bunching, an instability where the 

^ is important to have a model, which includes the kinetic 

■ steps move close to each other and form groups, called 

s_j ■ processes and the different effects existing on the crystal 

. 

surface. The classical model for the evolution of vicinal 



step bunches [HIIE]. 



surfaces was introduced by Burton, Cabrera and Frank The theoretical description of step bunching instabili- 

(BCF) [IJ. It is based on the observation that the kink ties within the framework of the BCF- model and its exten- 

sites are those positions at the surface steps where the sions has been the subject of much recent interest [HllTllHMni 

exchange between the adatom layer on the terraces and [T0l[m[T2lll3lfT4] . Here we focus specifically on the effect of 

the solid phase takes place. On the mesoscopic scale the non-conservative processes on the non-linear evolution of a 
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step train. As we reported in [14) for the problem of subli- 
mation in the presence of Ehrlich-Schwoebel (ES) barriers 
|15[|16j . non-conservative terms violating volume conser- 
vation in the co-moving frame arise generically from the 
interplay of sublimation and step-step interactions, and 
cause the interruption of the coarsening of the growing Fig- 1- (Color online) Sketch of the elementary processes in 




threshold [5]. In this setting the non-conserved dynamics 



bunches or splitting of a large bunch into several smaller the Burton-Cabrera-Frank model. 

bunches. In the present paper we expand this analysis to 

include the experimentally relevant effect of surface elec- . , , , , i , ^ • i • i 

IS described on large scales by the Benney equation, which 

displays either spatio-temporal chaos or an ordered array 

Latyshev and collaborators discovered that by changing r i , , • i?r7nrmn rr^, • 

of bunches, but no coarsening |281I29) . ihis macroscopic 

the direction of the direct heating current, a vicinal Si(lll) , , . . . ■ , , 

behavior is consistent with the complex mesoscopic step 

surface switches between bunching and debunching [T3. . , , • , 

dynamics revealed m the present work. 

Additionally, they observed several distinct temperature 

The paper is organized as follows. First, we sketch the 

regimes. In the so called regimes I and III 51 the bunch- 
derivation of the discrete step equations for the case of 

ing instability occurs only if the heating current is ap- 

attachment-detachment limited kinetics and present the 

plied in the down-step direction. On the other hand, for 

result of the linear stability analysis for very large wave 

the same direction in regime II debunching occurs, and 

lengths. Additionally, for comparison, we write down the 

bunching requires an up-step current. Here, we consider 

corresponding equations for the case of growth. We then 

the first temperature regime, where the temperature is low 

discuss the dicrete equations and their continuum limits 

enough in order to neglect step transparency (the motion 

for two special cases, where the kinetic asymmetry be- 



of adatoms across steps) ^ ^ . 

tween ascending and descending steps is caused solely by 

Interrupted coarsening of electromigration-induced step an ES-effect or by electromigration, respectively. Finally, 

bunches in the presence of sublimation was previously ^e show the results of numerical simulations of the dis- 

observed numerically by Sato and Uwaha [6^, however crete step equations for the case with electromigration. 

a detailed analysis of the phenomenon was not carried 

out due to the complexity of their model. Other stud- , , , , 

^ ^ 2 Model 

ies have approached the problem within the framework of 

weakly nonlinear amplitude equations, which can be sys- We consider an ascending one-dimensional step train with 

tematically derived by an expansion around the instability step edges located at positions Xi (FiglT]). The starting 



M. Ivanov and J. Krug: Non-conserved dynamics during electromigation-induced step bunching 3 

point for the derivation of the equations of motion for the must hold, where Q is the cross section of an atomic site at 
steps is the balance equation the step. The labels +/— refer to quantities correspond- 



dt 



dx^ dx 



ing to the lower/upper terrace of a step. The fluxes f± 

^+F = (1) 

T depend on both the difference of the adatom concentra- 



for the concentration of adatoms ni{x,t) on the j-th ter- 
race of width li = Xi-i-i — Xi. Here Dg is the surface diffu- 
sion constant, and r is the average life time of an adatom 
before desorption. Together these two quantities define 



the diffusion length Ijj ~ y/DsT, which sets the scale of 
spatial variation for the adatom concentration. The ter- 



tion n{x) compared to its equilibrium value n^.q and on 
the attachment/detachment to the steps with kinetic co- 
efficients k± . If the condition fc+ > fc_ is fulfilled we speak 
about a standard ES effect 15 . It induces an asymmetry 
in the concentration profiles (x) quantified by the asym- 
metry parameter 



race is losing adatoms due to desorption at rate 1/r and — ^+ _ ^- ^+ (o^ 

gaining adatoms because of deposition with constant flux 
F. The applied direct heating force F^i causes a drift of 
adatoms for which we can use the Einstein relation and 
define a further length scale, the so called electromigra- 
tion length /^^ = ksT/Fei [TS1I5]. As is common in the 
field, we assume that the adatom concentration adjusts 
instantaneously to the slowly moving steps, an assump- 
tion that is know as the quasistatic approximation and 

amounts to setting dtni(x, t) — in An approach that 

Aim fl^ \ 

goes beyond this approximation was recently presented by — — - —g ^ — ) — : gvi. (4) 



where l± = Ds/k± are called kinetic lengths. 

Apart from the attachment kinetics, a second effect 
incorporated into the boundary conditions ([2]) is the step- 
step repulsion. The equilibrium concentration rieq is de- 
termined by the chemical potential Afii at the ith step 
through the relation n^q « n'^qi^ + A^i/ksT), and Afj,i 
depends on the widths of the two neighboring terraces k 
and li^i according to pi[?T| 



Ranguelov and Stoyanov [241125] . 

The general solution ni(x) of the ordinary differential 
equation ([T]) can be specified using mass conservation at 
the steps as boundary conditions. A terrace of width I is 
bounded by two steps with positions x = ±1/2, at which 
the flux continuity conditions 



where I is the mean terrace spacing and 5 is a dimension- 
less measure for the strength of repulsion between the 
steps [1131]. 

3 Step equations of motion 



n 
n 



— fn 

dx 

dn 

a 

ox 



^ Using Eqs. pi2l4p we find the concentrations ni{x) for 

= +fc_[n - Tie,], X = --, 

^ all terraces. The velocity of the i—ih step is then given 

+ 1^ '^1' 2 ' by the superposition of the fiuxes coming from the two 
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neighboring terraces as dxi /dt = + f-^-. Since the non- 
conservative terms of primary interest here arise from sub- 
hmation, we discuss separately the hmiting cases of pure 
subUmation (F = 0, 7 > 0) and pure growth {F > 0, ^ ~ 
0); of course, in a typical experimental setup both pro- 
cesses may proceed simultaneously. For the case of pure 
sublimation and in the limit /^^ ^ Id, we obtain the 
non-linear system 



dt 



f{i--i+) 
2 



I-+1+ 







1 

' Id 


Ci-1 


1% 


1 /l.-l 




2 


-f 1 


Si-1 + 





(5) 



where R — mP^^QDs, Si = sinh(^i/Z£)), Ci = cos]i{li/l]j) and 
7i = 1 + gvi. The result ([5]) contains all four length scales 
and illustrates the complicated functional dependence for 
a simple one-dimensional step train. 

To simplify these expressions we use the approximation 
of attachment-detachment limited kinetics, l± ^ I 

[5^155] . After some calculations along the lines of [13] we 
arrive at 



r: 



dt ' [ 2 2 

- [in + It+i) h - in + li-i) h- 



(6) 



Here a second asymmetry parameter 
bei = - - 



1- + U 



(7) 



incorporating the strength of electromigration has been 
introduced. Re — {n'^gnDs)/ljj — n^gil/r is the constant 
rate with which the surface changes volume in a unit time 



(in the absence of non-linear, non-conservative terms, see 
below) and U = + 

4 Linear stability 

Equations similar to © can be derived when the sur- 
face is subject to a growth flux but sublimation is absent 
(F > 0, 7 = 0). In that case the factor 7^ in front of 
the square bracket on the right hand side of which 
depends nonlinearly on the step coordinates, is replaced 
by the constant 

7 = -5f- (8) 



Analogous to the problem considered in (14j . this implies 
qualitatively different instability conditions for growth and 
sublimation. Performing a standard linear stability analy- 
sis, in the limit of large wavelength perturbations we find 
the instability conditions 

6''"'' = 2bei +bES >Gg for subhmation (9) 



Ft 

= 2bei TrbES > for growth. 



(10) 



In the case of growth step bunching merely requires the 
compound asymmetry parameter b^^ to be positive, whereas 
for sublimation the corresponding quantity fe^"'' needs to 
exceed a positive threshold value 6g. This is an impor- 
tant consequence of the qualitatively different contribu- 
tions to the balance eq. (HJ that arise from desorption and 
deposition, respectively. Note that in a general situation 
the instability conditions (|9ll0p can be combined into the 
form 6 = (1 — Ft /■nPeq)bES + 26ei > 65, which was already 
obtained in |12) . 
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5 Conservative and nonCOnservative dynamics under summation over i, their contribution is conservative 

and the only non-conservative contributions in eq. (jlip are 



Beyond the linear stability properties, a fundamental dif- 
ference between the scenarios of pure growth and subli- 
mation is that the surface dynamics is conservative during 
growth but not during sublimation [51114). Here conserva- 



the first terms multiplied by ji. 

For comparison, setting b^i = eq.® reduces to the 
equations derived in [H] . 



tive dynamics implies that the rate of volume change of , dx, t, , , , 

the crystal, obtained by summing the equations of motion gbss 

^ ik - k-i) H ^ {i^ik~i - i^ik) ■ (12) 

over all steps Xi , is independent of the surface configura- 
tion [34]. Indeed, replacing the 7; in front of the square The difference between the two cases is that the terms 
brackets on the right hand side of © by the constant (© proportional to gKj on the RHS of eq. ^ do not cancel 
and summing over i, one readily obtains = -FHL, ^nder summation with respect to i, and thus are non- 
where L is the total length of the crystal. conservative. As will be shown in the next section, this 
It is instructive to compare the structure of the non- gj^es rise to distinct contributions in the continuum limit, 
conservative contributions induced during sublimation by 
the configuration-dependent factors ji in ([5]) for the two 

step bunching instabilities driven by electromigration and ^ Continuum equations 

by an ES-effect, respectively. First we neglect the ES- 

effect, setting bss ^ 0, which simplifies © into the form In previous work a systematic method for deriving con- 

j tinuum equations of motion from the discrete step dy- 

namics was developed [^[TU] which was applied to the 

-hel{k-k-l) (11) ^ 

model (jl2l) in |14 . Briefly, the method can be seen as 

+ 2 (^*^»-i v^li -\- v^-il^-i v^j^ili) . a kind of Lagrange transformation [35] which replaces the 

The second group of terms on the RHS of eq. ^ with 'Lagrangian' dynamics of particle-like steps by the 'Eu- 

prefactor U arises from equilibrium step-step interactions lerian' evolution of the step density m{x,t). The latter 

and stabilizes the regular step train, while the third group in turn defines a continuous height profile h{x, t) through 

of terms describes the effect of electromigration, which is "^(^^j — fi- 

stabilising or destabilising depending on the sign of b^i, Here we wish to compare the two instability mecha- 

i.e., the direction of Fe/. The last group of terms arises nisms described by eqs. dTTT) and ((HI), respectively, on the 

from the interplay of electromigration and step-step in- continuum level. Following the procedure outlined in [14j 

teractions. Since the terms in this group cancel pairwise for both models, we find that the continuum evolution 
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equation takes the general form 

3grn^ m' 3U (m^) 



dh d_ 
dt dx 

39 {m^y 



2 6m3 



2m 



-Jb 



1 = 



(13) 



where primes denote spatial derivatives. Here time t is 
rescaled by Re, length x by the average step distance I, 
and height h is measured in units of the monoatomic step 
height. The terms inside the square brackets on the LHS 
are conservative, and the non-conservative contributions 
are collected on the RHS of ea. ((T5| . The two models 
and (|12l) differ in the form of the contribution Ji, to the 
conserved surface flux, and of the non-conservative term 
<Pi,. Labeling the contributions due the ES-effect by ES 
and those due to electromigration by e/, respectively, the 
conservative terms are given by 



tES ^ES jel ^ei I „ , / 



and the non-conserved contributions are 



ES _ "ig^ES {m^ 



2m 



0. 



(15) 



As was discussed above, the terms in eq. ljlip proportional 
to gbgi give rise to a conservative contribution, whereas 
the terms in (I12p proportional to gbEs contribute to the 
non-conservative part of the continuum equation. 

In earlier work based on the continuum approach [9l 
110) the non-conservative contributions were generally ne- 
glected because of the smallness of g [13] , and it was there- 
fore concluded that step bunching phenomena induced by 
electromigration and by the ES-effect belong to the same 
universality class [51l36j . However, it has subsequently be- 



come clear that small non-conservative terms may qualita- 
tively change the nonlinear dynamics of surface steps [13] , 
and the fact that these terms are of different form for the 
two instability mechanisms implies that their equivalence 
needs to be reexamined. In the following we therefore ex- 
plore the nonlinear behavior of the electromigration model 
(jlip using numerical simulations. 

7 Nonlinear step dynamics 

Numerical simulations of eq. (jlip were carried out using 
an odeint-type procedure |37j for systems of M steps with 
periodic boundary conditions. We consider the following 
ranges for the four independent parameters of the model: 
bei e [0,0.5], U e [0,0.5], g G [0,0.1] and M < 100. 
Another degree of freedom is provided by the choice of 
the initial condition. In general, we start the simulations 
with two types of initial step train configurations: either 
a randomly disturbed equidistant step train, or an ini- 
tial shock of closely spaced steps and a single large ter- 
race. Step trajectores are shown in the co-moving coor- 
dinate system Xi{t) = Xi(t) — It, and we normalize both 
the height of the (monoatomic) steps and the average ter- 
race width / to unity. The time t is rescaled by Re and 
we measure the integration time in time units (t.u.). For 
the description of the bunch geometry we use two mea- 
sures: the maximal slope mmax = maxi{m.i} and the min- 
imal curvature = nrini{Ki}, where = and 
Hi = — 8(^i+i — li)/{li+i + k)^ respectively. A step is de- 
fined to belong to a bunch, if its distance to the next 
closest step of the bunch is smaller than ^ = 1 . 
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Fig. 2. (Color online) An example for the splitting of a large bunch in a system of 80 steps with parameters hei = 0.4, g = 
0.05, U — 0.05. a) Time evolution of one of the steps. The inset shows a blowup around the onset of splitting. The period of 
oscillations prior to the breakup reflect the time required for the step to once traverse the (single) bunch, b) Plot of all step 
trajectories between 11280 t.u. and 11470 t.u. c) Comparison of the profiles after 6000 t.u., 11500 t.u., 13000 t.u. and 15000 
t.u. d) Time evolution of the globally maximal slope, the globally minimal curvature and the number of bunches. 



An important consequence of the non-conservative char- 
acter of the dynamics is the phenomenon of anti-coarsening, 
where an initial large step bunch splits into smaller bunches 
[14] . In fig. [2] we show an example of this behavior for 
parameter values b^i = 0.4, g = 0.05, U = 0.05, and 
M = 80. Figure [5^) displays the movement of the first 
step in the train. After 11180 t.u. the large bunch splits 
for the first time into a large and a very small bunch, which 



is reflected in a clear change in the velocity of the step; 
smaller bunches move faster. Focusing on the splitting re- 
gion one finds that the initial small bunch disappears again 
after 4 t.u., an event which is repeated at 11258 t.u. with 
a small bunch life time of 19 t.u. In fig. [^b) we show the 
trajectories of all steps in the train in the time window 
11280- 11470 t.u. At 11333 t.u. the small bunch appears 
for the third time and at 11415 t.u. a third bunch arises 
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Fig. 3. (Color online) System with [7=0.2, M=40 and fluctuating initial condition, a) Stability/instability diagram, showing 
the number of bunches in the final state, for different combinations of b^i and g - □; 1 bunch, •: 1 or 2 bunches, x: 2 bunches, 
below the line b^i = 3g: stability, b) Time evolution of the maximal slope rrimax with bei— 0.35 and (from top to bottom) gi=0.00, 
0.01, 0.02, 0.05, and 0.09. 



for the first time. Tlie switcliing between different num- 
bers of bunches continues until the step train relaxes into 
four bunches, as seen in the height profiles in fig. [2J;). 
Figure [2ji) shows the corresponding time evolution of the 
minimal curvature and the maximal slope rrimax ^ 

along with the number of bunches. Both Umin and rrimax 
show clear changes in the region of the splitting; however, 
whereas the maximal slope remains essentially the same 
after the splitting event, there is a significant decrease in 
the minimal curvature. 

In fig. 13] we summarize results obtained from simula- 
tions starting from a randomly disturbed equidistant step 
array of 40 steps. Figure [3ti) shows the phase diagram 
of the system in the g-b^i plane at fixed U = 0.2. Be- 
low the line bd — 3g the system is linearly stable and 
rrimax — rn{x) = 1. Above this line we see three qualita- 
tively different types of long-time behavior: steady solu- 



tions with one single bunch, with two bunches, and time- 
dependent solutions that switch periodically between one 
and two bunches. In fig. [5]d) we plot the time evolution 
of the maximal slope rrimax at five points along the line 
bei—0.35. For g—0 we see the usual coarsening behavior, 
in which the number of bunches decreases in a step-wise 
fashion until a single bunch configuration is reached and 
the system relaxes to a stationary periodic state with a 
clearly bounded maximal slope. The remaining temporal 
periodicity of rrimax is due to the permanent step exchange 
between the front and back end of the bunch (see also inset 
in figl^h-))- Increasing g the maximal slope decreases (while 
still maintaining the single bunch configuration), until at 
g = 0.05 the regime of periodic switching is reached, lead- 
ing to a complex periodic pattern in rrimax- 

Finally, in fig. |4] we plot the behavior of the maximal 
slope as a function of the number of steps for two differ- 
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8 Conclusion 

In this work we have extended the non-conservative step 
bunching model presented in [T3] to include the effect of 
electromigration. The model applies to the first of the ex- 
perimentally observed temperature regimes on the Si(lll) 
surface, where step transparency can be neglected. The 
general step equations of motion incorporating sublima- 
tion, the Ehrlich-Schwoebel effect, electromigration and 
Fig. 4. (Color online) Dependence of the globally maximal step-step interactions were derived from the classical BCF 
slope m,nax on the number of steps M for g = [/ = 0.04, model in the quasistatic approximation. For the case of 
bei = 0.2 and for large and small noise amplitude in the random attachment-detachment limited kinetics we compared the 
initial condition. equations for growth and sublimation. In previous pub- 

lications non-conservative contributions were 

neglected, because of the experimentally small pref actor g 
Plll4j . Those terms were now taken into account and some 
important consequences were identified. First, on the level 
ent amplitudes of the initial disturbance, (7=[/=0.04, and of linear stability analysis, they shift the instability condi- 
bei— 0.2. Here rrimax is the global maximal slope, measured tion on the dimensionless asymmetry parameter b by 6g, 
as its largest value for the last 500 t.u. of the simulation, as was first pointed out in [T2]. This shift is present in the 
We see that m,„cia; generally increases with M, but this be- case of sublimation, but not in the case of growth [T3]. 
havior is interrupted by downward jumps every time the Moreover, in the case of sublimation the structure of the 
number of bunches that can fit into the system increases non-conservative terms differs depending on the underly- 
by one. This shows that the existence of stationary solu- ing mechanism inducing the asymmetry between ascend- 
tions with multiple bunches can be seen as a consequence ing and descending steps. This leads to different contin- 
of the fact that, in the presence of non-conservative pro- uum equations for step bunching caused by an ES-effect 
cesses, the maximal slope is bounded from above (l4j . Near or by electromigration, respectively. 

the transition between different numbers of bunches the Nevertheless, the numerical integration of the discrete 
system behavior depends very sensitively on the ampli- step equations for the case with sublimation and electro- 
tude of fiuctuations in the initial configuration, an effect migration reproduces qualitatively the results of |14) . The 



that is particularly pronounced around M = 70. 



non-linear, non-conservative terms supply a richness of 
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